Conformation of Circular DNA in 2 Dimensions 
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The conformation of circular DNA molecules of various lengths adsorbed in a 2D conformation 
on a mica surface is studied. The results confirm the conjecture that the critical exponent v is 
topologically invariant and equal to the SAW value (in the present case v = 3/4), and that the 
topology and dimensionality of the system strongly influences the cross-over between the rigid 
regime and the self-avoiding regime at a scale L w 8 i p . Additionally, the bond correlation function 
scales with the molecular length L as predicted. For molecular lengths L < 5 £ p , circular DNA 
behaves like a stiff molecule with approximately elliptic shape. 



From a theoretical and experimental perspective, the 
physics of linear polymers is well-established llj and its 
predictions were also confirmed for linear DNA [1, H, 0, HI • 
On the contrary, topologically constrained polymers, like 
circular or knotted DNA, are still challenging, and the 
topological effects on the static and dynam ic p roperties 
are a matter of debate @, 0, S, H EE EH, G3- A con- 
vergence has been reached on the theoretical value of the 
critical exponent v describing the divergence of the ra- 
dius of gyration with the length of circular and knotted 
polymers @, EH, E3] , and fixes this value to the SAW ex- 
ponent. Recently, experimental evidence has been found 
that the scaling exponent v for knotted DNA in 3D cor- 
responds to the 3D linear SAW value fl5j | . For circular 
DNA in 3D, diffusion measurements have pointed in the 
same direction [T3|. Real polymers, like DNA, have a 
mechanical stiffness which strongly influences the scaling 
properties as well as the bond correlation function, mak- 
ing the theoretical treatment more complex. Semiflexible 
polymer models wereput forward in order to include this 
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into the theory [171 . Il8t [19| , and showed the existence of 
a gaussian regime between the rigid and the self-avoiding 
regimes Many theoretical investigations concerning 
the properties of circular molecules were motivated by the 
observation of DNA in closed form 120, 21, 



23. 24 



These studies are relevant for biology because DNA fre- 
quently adopts circular form in vivo [25J. Furthermore, 
the large persistence length l v of DNA permits to test 
semi-flexible polymer theories and the effect of correla- 
tion along the molecule. 

In this Letter, we present the first results for circular 
DNA in 2D spanning lengths from L/£ p = 2 up to 120, 
where L is the plasmid length. We show that the gaus- 
sian regime is suppressed in 2D in favor of a direct transi- 
tion from the rigid rod behavior to the self-avoiding walk 
behavior indicating a strong effect of the dimensionality 
and of the topology on this transition Q . The scaling of 
(Rq) is confirmed, and the internal end-to-end distance 
(£ 2 ) is well-described by the Bloomfield-Zimm formula 
[2 11 ] . The correlation function scales as predicted by the- 
oretical studies 22, 23| with the critical exponent for a 
SAW in 2D. For the short plasmids, L < 5 l v a rigid rod 



FIG. 1: AFM images of DNA plasmids deposited on Mg 2+ - 
mica. (a) minicircle-1 and -2 (350 nm scale) (b) pUC19 (2 
/jm scale), (c) pBR322 (3 (jm scale), (d) pSHl monomeric, 
dimeric, and trimeric forms (2.5 /jm scale). 



behavior for the tested quantities was observed. 

Plasmids pUC19 (2686 bp), pBR322 (4361 bp), and 
pSHl (5930 bp) were nicked to produce circular molecules 
without superhelicity. Plasmid pCU19 was treated with 
restriction enzyme Rsal to produce 3 different linear frag- 
ments, and using T4 DNA ligase mini-circles of differ- 
ent length were obtained. All DNA preparations were 
diluted in 1 mM Tris-HCl buffer, pH 7.8, to a final 
DNA concentration of 1 /xg/ml. MgCl2 was added in 
the DNA solution to a final concentration of 5 mM. A 
10 /il aliquot of the DNA solution was deposited onto 
freshly cleaved mica, and incubated for 10 minutes at 
room temperature. The sample was then rinsed with 
nanopure (Ultra High Quality) water (USF Elga, High 
Wycombe, England) and dried with air. Images were 
collected using a Nanoscope Ilia (Veeco Inc, USA) op- 
erated in tapping mode in air. Ultrasharp, non-contact 
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silicon cantilevers (NT-MDT Co., Zelenograd, Moscow, 
Russia) with a nominal tip radius of < 10 nm were driven 
at oscillation frequencies in the range of 150 to 300 kHz. 
During imaging, the surface was scanned at a rate of one 
line per second. Images were simply flattened using the 
Nanoscope III software, and no further image process- 
ing was done. The contour of the DNA molecules was 
extracted with Ellipse [26j |. 

From the AFM images shown in Fig. [TJ it is readily 
clear that the molecules exhibit no self-crossings, indicat- 
ing that they equilibrate in a 2D conformation. This find- 
ing is supported by an analysis of the persistence length, 
£ p , using the bond-correlation function for semi-flexible 
polymers in 2D: 



(cos 9(s)) = exp(~s/2£ p ) 



(1) 



where 9 is the angle between the tangent vectors to the 
chain at two points separated by the distance s. A proper 
fitting can only be done for molecules with L > 10 £ p 
[3pj ] . below this length, the mechanical stiffness imposes a 
disk shape. Fitting the experimental correlation function 
with Eq. (JXJ) for small separations s, values of £ p close to 
50 nm for pUC19, pBR322 and pSHl (see Table [Q) were 
found, in good agreement with the one for linear DNA Q , 
indicating that the molecules do adopt an equilibrated 2D 
conformation. Using the 3D equivalent to Eq. |lj one 
would overestimate £ p by a factor of 2 0, [2t| • 

In Fig. 03 the bond correlation functions for the entire 
range of segment separations s are depicted for plasmid 
lengths L > 10 l v . Unlike the case of linear chains, after 
the initial exponential decay, circularity imposes negative 
values with a symmetrical minimum at s — L/2. In the 
case of a long flexible chain, the bond correlation func- 
tion was studied in several theoretical papers [22j,l23( . For 
a random walk, the correlation would be uniformly dis- 
tributed over all segments [22j. However, self-avoidance 
alters this distribution, and farther spaced segments will, 
on average, have a more negative correlation depending 
on the total length. A scaling relation for the bond cor- 
relation function was suggested by Baumgartner [23j; 



-2u-2 



(2) 



where is a scaling function depending on the relative 
length s/L, v the excluded volume, and v the critical ex- 
ponent pertaining to the scaling of the radius of gyration. 
Rescaling our data according to Eq. ([2]), and using the 
critical exponent for a SAW in 2D [y = 0.75), the data 
collapse on one universal curve as shown in Fig. OA. The 
fact that v ~ 0.75 results in a good data collapse is a 
further confirmation of the 2D equilibration of the DNA 
plasmids on long length scales. 

For the two smallest molecules (L — 1.6 £ p and 4.6 £ p ), 
the correlation function (depicted in Fig. [3]) is closer 
to that of an ellipse of eccentricity 0.67 [10|. A small 
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FIG. 2: (a) Bond correlation function scaled according to Eq. 
p| with v — 0.75. (b) Unsealed bond correlation function. 
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FIG. 3: Bond correlations for the smallest molecules. The 
theoretical correlation function for a circle (dashed) and for 
an ellipse (continous) with eccentricity 0.67 are also drawn. 



deviation around s = L/2 still remains, whose origin is 



not yet clear 28| 

Next, we determined the radius of gyration (Rq) as 
a function of the molecular length L in order to test 
the scaling relation (Rq) ~ L 2v . To ensure good statis- 
tics, the average for each molecular length was taken over 
about 100 molecules (except the case of pSHl-trimer with 
only 25 molecules) . In Fig. |4] we show a log-log plot of 
(R G ) vs. L with power law fits of the data, giving directly 
the value of the critical exponent v. Two distinct regimes 
appear on this graph. A fit through the two data points 
for smallest L gives a value of v = 0.86 ± 0.03, expressing 
the fact that for these almost circular molecules the ra- 
dius of gyration Rq is equal to the radius of a circle, thus 
Rg = L/2ir. In Fig. [l]the dashed line represents the ra- 
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FIG. 4: Radius of gyration as a function of the plasmid length. 
Fitting with (R%) ~ L 2v gives j/ = 0.86 ± 0.03 for the small 
scale regime and v = 0.75 ± 0.04 for the large scale regime. 



TABLE I: Properties of the investigated 
length L is calculated using 0.34 nm per 
are in nanometers. 



DNA molecules. The 
basepair. All lengths 



DNA type 


(bp) 


L 


L/l p 


£ P 


Let 


V 


mini fjl 


241 


82 


1.6 


N.A. 


83±1 


0.90±0.01 


mini J 2 


676 


230 


4.6 


N.A. 


241±1 


0.91±0.01 


pUC19 


2686 


913 


18.3 


52.5±5.1 


937±13 


0.85±0.07 


pBR322 


4361 


1483 


30 


49.1±4.6 


1543±5 


0.82±0.01 


pSHl 


5930 


2016 


40 


53.5±4.2 


2110±10 


0.81±0.01 


2 x pSHl 


11860 


4032 


80 


49.1±2.8 


4125±12 


0.76±0.01 


3 x pSHl 


17790 


6048 


120 


49.6±2.5 


6293±18 


0.75±0.01 



dius of gyration \Rq) for a perfect circle. For the larger 
molecules, we obtain a value v = 0.75±0.04, in very good 
agreement with the theoretical value for a SAW in 2D. 
This clearly shows that circular molecules exhibit scaling 
with molecular length, and that the scaling exponent is 
topologically invariant. 

Fig. [4] also suggests that there is no gaussian regime 
[y = 1/2) between the stiff and the SAW regimes as it 
was found for linear polymers [3| • This behavior could be 
due to the topology and the dimensionality of our DNA. 
In fact, Camacho et al. [|| described the scaling prop- 
erties of semi-flexible circular polymers in 2D, and did 
not find an intermediate gaussian regime. The reasons 
are that collisions between two strands are more likely in 
2D, and that circularity imposes to the polymer strand to 
come back on itself: both effects increase the self-avoiding 
interaction on short distances. From Fig. |4]the cross-over 
length (= the length at which the self-avoidance mani- 
fests in the radius of gyration) was estimated by locating 
the crossing of the lines with slope v = 1 and v = 3/4. 
This gave L = 400 nm ~ 8£ p , agreeing well with the 
results in Ref. tgl. 
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FIG. 5: The internal end-to-end distance yj (£ 2 ) as a function 
of contour length s for seven lengths of DNA plasmids. Fits 
are according to Eq. ((3}, and the results are summarized in 
Table U 



An alternative method to determine the critical expo- 
nent v uses the separations of segments within individ- 
ual molecules. Although this quantity was measured for 
linear chains hardly any experimental data exist for 
ring polymers. Numerical values for v were determined 
by Ref. 14(. For a linear molecule the scaling of intra- 
chain separations equals that of the end-to-end distance, 
and thus a larger contour length implies a larger sepa- 
ration. However, for a circular molecule the topological 
constraint limits the intrachain separations to L/2, and a 
specific expression must be used. For a ring polymer with 
excluded volume effects, the internal end-to-end distance 
(£ 2 (s)) is given by 



<£ 2 (s)> 



s 2 v (L-sf 
+ (L- s) 



.2;/ 



2v ' 



(3) 



where s the contour length. The data for the intrachain 
separations are shown in Fig. [5] Several interesting ob- 
servations can be made. As the molecules are circular, 
the value of (£ 2 ) increases only up to contour lengths 
s ~ L/2. For larger values, the end-to-end distance de- 
creases back to zero. In the case of the largest molecules, 
for contour length below t p , a line with slope one de- 
scribes the data well, the reason being that the molecules 
on this length scale are rod-like. Increasing the contour 
length beyond t p , the self-avoiding regime is entered. 
Here, the molecules are highly flexible, and are well- 
described by SAW walk with circular constraint. Fitting 
this large-intrachain separation regime with Eq. ([3]), we 
obtain values for the critical exponent i/, and the total 
contour length L as shown in Table HI For the longest 
molecules good agreement with the expected molecular 
length as well as the exponent is found. However, for 
short molecular lengths the fitted value of v deviates sig- 
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nificantly from 0.75, while L still is in good accord with 
the expected length. As the relative importance of £ p to 
L increases for small chains, the available size of the SAW 
regime becomes limited and the internal end-to-end dis- 
tance £ is better described by the one corresponding to a 
circle, namely £ 2 (s) = 2(L/2n) 2 (l - cos(27rs/ L)) or with 
a similar expression for an ellipse. For the two smallest 
minicircles, only the stiff regime is visible, because the 
molecules are too short to enter the SAW regime and 
thus the critical exponent has a value of v ~ 0.9. 

The internal end-to-end distances (£ 2 (s)) is well-fitted 
with Eq. [3J and therefore should follow a universal curve 
if V (£ 2 ( s )} divided by V and the contour length s by L. 
Indeed this is the case and the rescaled plot is presented 
in the supplementary material [30j ]. 

In conclusion, we have shown that circular DNA 
molecules adsorbed on Mg 2+ -mica equilibrate in 2D con- 
formation, and we were able to confirm several theoretical 
predictions concerning ring polymers. The bond corre- 
lation function obeys well to the predicted scaling [23| 
within reasonable error when long plasmids are consid- 
ered L >> l p : this scaling seems to be reached only for 
the two longest plasmids (L = 80 & 120 £ p ). In partic- 
ular the bond correlation function (not presented here, 
but in Ref.[30|]) decays exponentially for contour lengths 
s < £ p , and this leads to the determination of the per- 
sistence lengths reported in Table fl] Beyond the initial 
exponential part, nothing is known about the exact de- 
pendence of correlation function upon s for s > £ p : this 
dependence is not trivial as it can be seen in Fig. [2] 
It must be noticed that the scaled correlation functions 
cross each other at the reduced length s/L = 1/tt. Addi- 
tional theoretical studies are needed to understand this 
behavior, especially in the case of polyelectrolytes like 
DNA. For short plasmids (L < 5 £ p ), the correlation 
function can be better described by that of an ellipse of 
eccentricity e = 0.67, whose origin could be fluctuations 
[h| or flexibility [H, [3l|] . This fact rises questions about 
the reasons of the deviations from the circular shape and 
should further investigated both experimentally and the- 
oretically 1(J 2i|. Furthermore, from (Rq(s)) vs. L, a 
cross-over from stiff to flexible SAW behavior is observed 
without an intermediate gaussian regime (Fig. QJ. This 
transition takes place at L ~ 8 £ p , a number that seems 
due to the topological constraint of circularity Q. The 
critical exponent v was found to be v = 0.75, and ex- 
hibits topological invariance. Similarly, from the internal 
end-to-end distance (£ 2 (s)) vs. s (Fig. [3]) we confirm 
the predictions by Ref. [2l| . Interestingly Eq. [3] is valid 
for stiff as well as flexible circular molecules. In sum- 
mary, the available plasmid lengths permitted to cover a 
large range of conformations, from totally stiff chains to 
SAW chains. The importance of stiffness was shown on 
molecules with a length of a few £ p . 
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